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§ 1. Introduction, 

In a previous paper"^ the authors, with others, have described observations of the 
angular oscillations of the axis of a 3-inch shell over the first 600 feet from the muzzle 
of the gun, and from an analysis of the observations have obtained information about 
the forces due to the air. In the experiments, shells were fired from two guns giving 
different degrees of axial spin to the shell. While the shells fired from the gun giving 
the more rapid spin were all stable, most of the shells from the other gun were slightly 
unstable, this condition being shown by the much larger maximum yawf developed. 
These unstable rounds were not analysed in (A) as no suitable method of doing so had 
then been devised. | The analysis of these rounds, about one-third of the number 
fired, forms the subject of the present paper ; the results confirm those of (A) and 
provide some additional information. 

The information as to the force system obtained from the stable rounds was 
confined to yaws up to 7 degrees or perhaps 10 degrees ; by analysis of the unstable 
rounds this information is extended, though in a fragmentary manner, over the region 
up to 35 degrees of yaw. On the other hand, no information has been derived from 
the observed damping of the unstable rounds. The observations are, in respect of 
the damping, clearly in qualitative agreement with the theory and results of (A), but 
no method has been devised of making a quantitative analysis of the damping. 

The force system on a model shell was also determined at low velocity in the wind 
channels of the National Physical Laboratory. The results are shown in fig. 2 of (A) 

* ''The Aerodynamics of a Spinning Shell," 'Phil. Trans.,' A, vol. 221, p. 295 (1920). This paper 
will be cited as (A). The experiments analysed here and in (A) were carried out for the Ordnance 
Committee, and the results are published with their sanction. 

f The '' yaw " is the angle between the axis of the shell and the direction of motion of its centre of 
gravity. 

{ As will appear later, the ordinary solution in elliptic functions of the equations of motion of a top 
is not adequate for this purpose in the case of large yaws. 
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and have been used in both, papers for extending the results down to low velocities, 
as in figs. 1 and 2, here. 

Shells of four different types, I.-IV. were used. Types I.-III. were of the same 
external shape (form A), with three diiferent positions of the centre of gravity. 
Type IV. was of a different external shape (form B). The details are given in (A).* 

The experimental data, which have already been discussed in (A), consist of the 
mass, principal m.oments of inertia, and position of the centre of gravity for each type 
of shell ; rough values of the forward velocity over the whole range of 600 feet from 
the muzzle of the gun ; the spin of the shell, deduced from the rifling of the bore ; the 
yaw and orientation of the shell's axis at a number of points along the range, deduced 
from the shape and orientation of the holes punched by the shell in cardboard targets. 
These cards were set up at intervals of about 60 feet for all the unstable rounds, and 
it appears from figs. 12 of (A) and figs. 3 and 4 of this paper that they were close enough 
together for satisfactory curves to be drawn through the observed points representing 
the variation of the yaw S and its azimuth <^ over the whole range. 



§ 2. The Equations of Motion. 

It is convenient to recapitulate the notation of (A). Suppose that OA denotes the 
direction of the axis, OP the direction of motion of the centre of gravity of the shell ; 
then AOP = (5, and ^ is the angle that the plane AOP makes with a fixed plane 
through OP. M (= M sin ^) is the couple in the plane AOP which tends to increase (5, 
A, B and N are the principal moments of inertia and the axial spin of the shell, and we 
write Q = AN/B. Then the equations of motion will be taken in the formf 

r + ^'' sin' ^+ f % d cos (5 - E, . . (1) 

0' sin' (^ + 12 cos (5 = P, (2) 

where E and F are constants. The various assumptions underlying equations (1) 
and (2) are discussed in detail in (A). If ^x is constant the equations are, of course, 
of the same form as the ordinary integrals of energy and angular momentum for a 
spinning top, and the complete solution in elliptic functions is standard. 

When M is an arbitrary (odd) function of § the top solution no longer applies, but 
a solution in elliptic functions is still possible if M has the form X sin (^ {1— Y (l^cos 8)}, 
where X and Y are constants. This more general form allows the first two terms in 
the expansion of an arbitrary M to be catered for and can represent M adequately 

* Log. cit., p. 316 and fig. 6. See also fig. 1, here. 

I (A), loc, cit., p. 334, equations 3.404, 3.405. For tlie underlying assumptions see (A) Part I., pp, 301 sqq., 
311 sqq. These equations are, strictly speaking, not referred to fixed axes, but are approximate equations 
referred to axes changing direction with OP. Dashes denote difierentiations with respect to the time t. 



THE AERODYNAMICS OF A SPINNING SHELL. 229 

over a wider range of values of §. By suitably adjusting X and Y, which define the 
couple, and the initial conditions, a curve showing the variation of § with the time 
can be obtained which agrees closely with observations over a complete half -period, 
so that the above expression for M appears to be adequate up to values of § of 35 degrees.^ 
Our original approximation with Y == fails in general when (^>10 degrees. The 
observed curves suffice to determine X and Y for each round, and as observations 
w^ere taken for a number of different values of the muzzle velocity, M is determined 
by the experiments over a limited range, as a function of the two variables, v, the 
velocity of the shell, and §. 

In solving the equations of motion it is convenient to express the couple by means 
of the non-dimensional coejfficients s and q defined by the equation 

^^Ba^Alll{i^^q,{l^^o^§)}, . (3) 

It will appear that the motion with S permanently zero is stable or unstable 
according as 5>1 or 5<L For the rounds here analysed, s lies between 1-06 
and 0-83. 

In expressing the results in a standard form it is convenient to use a different 
non-dimensional coefficient /m, which is independent of the mass, moments of inertia, 
size and velocity of the shell, and depends only on the shape of the shell and the 
non-dimensional variables v/a and §, where a is the velocity of sound. This is defined 

by the equationf 

M^ pvV^mSfy,{v/a,§), (4) 

where p is the air-density, r the radius of the shell, and the quantities involved are 
expressed in consistent units. 

According to (3), /m is practically constant so long as (5 < 7 degrees, and the value 
of /m ('^/<^5 0) is strictly comparable with similar values obtained in (A) by analysis 
of the stable rounds on the assumption that/^ is independent of d 

§ 3. Final Results of the Experiment. 

Fig. 1 shows curves of /m {v/a, 0) as a function of v/a for the four types of shell 
corrected for the effect of the cards J. They are reproduced without alteration from 
figs. 4 and 5 of (A) and represent the results for the stable rounds. § The values 
derived from the present analysis of the unstable rounds are plotted for comparison ; 

* When tlie yaw exceeds 30 degrees the fit is less satisfactory {e.g., III., 11-13). 
t Log. cit., p. 302, equation 1.103. 
J § 10 below. 

§ The curve for type II. is not actually given, but the data for drawing it can be found in (A) (iig. 13, 
p. 352). 

2 K 2 
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they show remarkably good agreement between the new results and the old. This 
confirms the substantial truth of the whole theory ; in particular, the agreement of 
the results for the two twists of rifling verifies that the couple M is unaffected by a 
change in axial spin in the ratio 3 to 4. 




Fig. 1 . Values of the couple coefficient /m (vja, 0) and the normal force coefficient fj^ (v/a, 0) plotted 
against v/a, corrected for the effect of the cards. The full curves for/ai for types I. -IV. are reproduced 
without alteration from (A). The plotted points show the means of observations here analysed for the 
first and second half-periods, and the numbers against them show the number of rounds contributing to 
each mean. The group marked * was analysed in (A) by a different method with identical results. The 
origin of /^ for type IV. is displaced downwards 4 units. 

The full curve (V) for f^ (shells of form A) represents the complete results of the experiment. The 
dotted curve reproduces the partial results of (A). 

Curve l.—fu for 3-inch shells of form A with centre of gravity 4-20 inches from the base. (Type III.) 

Curve II. — The same Avith centre of gravity 4-73 inches from the base. (Type I.) 

Curve III. — The same with centre of gravity 5-08 inches from the base. (Type II.) 

Curve IV. — The same for form B with centre of gravity 4-965 inches from the base. (Type IV.) 

Curves of /^ sin §, the complete moment coefficient, considered as a function of 
both variables v/a and S corrected for the effect of the cards, are plotted against 
S in fig. 2. The information is somewhat fragmentary : in addition to wind channel 
results, values of /m, when ^> 10 degrees, are available for shells of type I -III. for 
two velocities near v = a. For type IV. (pointed shells) values of /m for large § are 
wanting in this region, but exist for two high velocities, and one less than a. 
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Some general information can be deduced from these curves. At high velocities 
(v/a = 2-0) and large values of § the couple coefficient actually falls below its low 
n 




5 10 15 20 25 30 35 40 

Fig. 2. Curves of tlie complete couple coefficient /m (vja, 5) sin 8 against S for various values of v/a 

(as shown against the curves). The curves stop at the greatest value of S for which observations are 

available. 

The origins of all the curves except for type III. are displaced upwards and can be recovered from the 

fact that all curves pass through their respective origins. The scale of S is in degrees. 

velocity (wind channel) value (fig. 2) ; the curve of the couple coefficient against yaw 
has here a large . curvature downwards. It is, on the other hand, almost straight in 
the region of the velocity of sound. 
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We may notice also the peculiar behaviour of rounds III. (11-13), fig. 3,"^ in which 
the motion seems to become quite irregular from round to round after the first 
maximum yaw is attained. As these are the only rounds in which a yaw of 35 degrees 
or more is developed, it seems likely that some peculiar change in the type of the 
airflow occurs at or about this angle, analogous to the phenomenon of " burbling. "f 

Eounds IV. (10-12) have the smallest velocity of any group fired ; the agreement 
of /m for these rounds with wind channel results for the whole range of §, figs. 1 and 2, 
is very satisfactory. 

The new values of /^ {v/a, 0) can be used to correct slightly the mean curves of /m, 
and from the modified curves the value of /n, the coefficient of the force normal to 
the shell's axis at small yaw, may be re-determined by the method used in (A), § 1-13. 
The re-determined curve for f^ is shown in fig. 1 . The values differ only slightly from 
the former values, and the main features of the f^ curve are fully confirmed. 

The final results of the complete experiment, stable and unstable rounds alike, 
corrected for the effect of the cards, have been combined together to give the final 
values of /^ {v/a, 0) for shells of types I. and IV. s^ndf^ {vja, 0) for type I. shown in 
Table I., which replaces the corresponding Table I. of (A)» 



Table I. — Final values of /^ {vja, 0) and f^ (v/a^ 0) for shells of type I. and/^ (v/a, 0) 
for shells of type IV., embodying the results of the whole experiment. The effect of 
the cards has been corrected for as far as possible, and this table supersedes Table I. 
of (A). 





v/a. 


1 

Shells of type I. 


Type IV. 












/m {v/a, 0). 


fi, {via, 0). 
3-34 


" [ 

8-95 




Wind cliannel 


8-57 




0-7 


8-6 


3-35 


9-0 




0-8 


9-05 


4-0 


9-7 




0-9 


10-35 


5-2 


11-1 




1-0 


11-55 


5-25 


11-75 




1-1 


11-55 


4-7 


11-6 




1-2 


11-25 


3-9 


11-4 




1-3 


10-9 


3-7 


11-25 




1-4 


10-55 


3-85 


11-1 




1-5 


10-3 


4-0 


11-0 




1-6 


10-05 


4-15 


10-95 1 




1-7 


9-85. 


4-3 1 


10-9 ! 




1-8 


9-65 4-5 


10-8 




1-9 i 


9-4 


10-75 1 




2-0 


9-15 







* Also (A), p. 349, fig. 12b. 

t This term is commonly applied to the sudden increase of turbulence behind an aerofoil at the critical 
angle. 
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^ A: The Solution of the EquafAons of Motion, 

We shall now solve the equations of motion (1) and (2), assuming that M is of form (3). 
We take first the case of rosette motion, in which zero values of § can occur, so that we 
may assume the initial conditions 

Eliminating ^' and writing sin ^S =^ y we get 

Ay'' = Q^{b^l-f)-,/+^{l-y^){l/s-iqf)} (5) 

The right-hand side is a cubic in ^^ whose roots are such that it may be written in 

the form. 

iy" = Aqif{h'+f){a'-f){r-i/'), (l//^ < 1< iK). ... (6) 

Formulae connecting* h, a and / with b, q and s may be obtained most conveniently 
by putting y^ =: i^ y^ =: o, and by differentiating with respect to y^ and putting 
y^ z=z 0. The resulting formulaa are 

4g(l+A^)(l-a^)(r-l)^l, . (7) 

. iqa'KV == b\ . (8) 

Aqia'f-ct'h'-h'D = ~b'+l/s~l (9) 

A solution of (6) is obtained by assuming, in the usual notation of Jacob ian elliptic 
functions, 

2 a^ Cli^ U r^^^ 

y z=— ~- (10) 

where the constant g and the modulus Jc of the elliptic functions remain to be determined. 
If we solve (10) for cn^ it and differentiate, we get 

— sn u en nan ii u — -_«A~^-i^-^^^--, 

{a'-gYY 
leading to 

Comparing this with (6) we may write u^ '=^ ± XO, X constant, and obtain 

f=^ay\ (11) 

^^'^'4^' • • ^^^^ 

^~ a^{i-g') ■ • ■ ■ ■ • ^ 

'•^ This a lias, of course, no connection with the velocity of sound. 
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Since y = when 1 = 0, we may write u = K.—\Qt, where K is the complete 
elliptic integral of the first kind to modulus k. Treating Qt as independent variable, 
the final form of our solution contains four constants a, \y k and g, which must be 
completely determinable in terms of s, q and 6, so that there must be an independent 
relation between a, X, k and g by which any one can be found when the other three 
are known. To determine them it would be necessary to solve (5), the original cubic 
in y^. To analyse the experiments, however, we have to solve the inverse problem 
of determining s, q and b when a, A, k and g are known. In practice g is small so 
that, as a first approximation, we may use the following simplified form of (10) : — 

^ = acn (K— XQi^), (14) 

where a, k and X may be treated as independent. By fitting a curve of this type to 
the curve of observed values of y (sin |(!) against 0^, we can determine the constants 
a, X and k. It is at once clear that a =: sin |a, where a is the maximum yaw, but 
we shall continue to call this constant a for shortness. The value of g can then be 
obtained from the identical relation in terms of a, k, X, and the curve re-calculated by 
formula (10) if gr is large enough to make it worth while to do so. After that the 
values of a, X and k could be re- adjusted and the process repeated. Theoretically, 
we could presumably arrive at the precise solution in this manner by a limiting 
process. Practically, in nearly every case, the first approximation with ^ = is all 
that is required. 

The values of s, q and b are given by simple formulae in terms of a, X, k and g. Prom 
(8) and (9) we get 

on substituting for q, W and/^ from (11)™(13), this reduces to 

-, /^ -I ..2/ o , cos^c(a^— 3g^)| t.^x 

1/^ — 1 = 4x -^ — cos 2/c-f- ^ .. K » ..... (Id) 

where k = sin k.^ In practice either a or cos k, or both, are small and g is of the same 
order as a ; the second term inside the bracket may then be neglected in determining s, 
in which case the value of g is not required. This, as we shall see, is really a consequence 
of the smallness of 6, its mean value in practice being about 0-015. For q write 
equation (7) in the form 

and substitute for q, ¥ and p in the right-hand side, getting 

4g = -_L^+^(sinV+22f£K=V}l (16) 

* This K will not be confused with the k of (A), loc, cit, p, 328, which is the damping coefficient depending 
on the cross-wind force. 



THE AEKODYNAMICS OF A SPINNING SHELL. 235 



Finally equation (8) gives 



¥ = !^!^!^ (17) 



The equation for g is obtained by substituting for q, h^ and/^ in (7), which, becomes 

4X'{l-^a'){a'-g'){a'-g'+{l-a;')F}-a'{l-g') = 0; . . . (18) 
this is a quadratic for g^^ whose solution may be written 

(/•^ = sin^ ia (l - *-^) , (sinia = a), (19) 

where 6 is given by 

cot 6 = XF cot' ia- (tan' ia)/4X. ....... (20) 

The ambiguity is settled in practice by the fact that g^ must be small if equation (14) 
is to be taken as a first approximation to the solution. In practice, as we have said, 
b is small. Valuable information as to the nature of the solution is, therefore, obtainable 
by considering its limit as & -> 0. This gives us a guide as to the actual relative order 
of all terms. 

Let us suppose then that h -> 0, that s and q are definite constants, and let us assume 
that g is of the same order as a, which by (19) and (20) must be the case unless X -> 0. 
Then equation (17) shows that X^a^ cos^ KJ{l—g^) ->0, and, therefore, in the limit, 

1/6>~1 = -4X'cos 2/c (21) 

Equation (21) shows that, ii s r^ 1, X->0 is impossible. Hence in alh cases {s ^ 1) 
our assumption as to g is justified and a cos /c^O. This also justifies our previous 
statement concerning (15). There are now two cases according as ^<1 or 5>1. 

Case (i). — s < 1. We are supposing that s—l is fixed, so that as & -^ 0, s—l is large 
compared to &. To satisfy the signs of (21) we must have /c > 45 degrees. This 
implies sin^ ^>\-> ^^ that (16) becomes in the limit 

, 1 , 4X' sin' K 

It follows that a -> is impossible, and therefore cos /c ^ 0, /c -> 90 degrees, and a and X 
tend to definite non-zero limits. The limiting forms for 1 js, q and g^ are easily found 

to be 

1/6^-1 = 4X^ (22) 

4g = — sec'|-a + 4X'cosec'|a, . . . ... . . (23) 

^'-sin'|a{l-(tan'|a)/4X'}, , . (24) 

VOL. COXXII, — A. 2 L 
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in which we have replaced a by sin \a. These formulae are good approximations in 

practice, when s is not too close to unity. The half -period QT"^ (— K/x) tends to 

infinity, but so slowly that no difficulty occurs in practice. Since a, the maximum 

yaw, does not tend to zero with &, the initial disturbance, we have here what may be 

called the unstable case. 

Case (ii). — 5>1. We now must have /<:<46 degrees and, therefore, as cos^ ^^h 

a->0 at the same rate as b ; this is the stable case in the usual sense. Equation (16) 

shows further that /c->0, and, therefore, by (15), \ has a definite limit determined 

by 

1-1/5 = 4X'. 

For given q equations (16)-(20) determine the limiting ratios of b : sin |a : k : g. 

The case in which s = 1 and 6 -> can be treated in a similar way. It is found that 
/c ^ 45 degrees and a and X both tend to zero like ^b. 

§ 5. Rounds ivith a Non-Zero Minimum Yaw, 

We shall only consider cases in which the minimum yaw /3 is small, and shall take 

as initial conditions 

' § = /3, §' == 0, <p' sin § - 61O. 

The equations of motion become 

0^sin^(5-Q6isin/3 + fi(cos/3--cos(5) = 0, (25) 



r' + cp''Bm'S--Q%'+f^dcos§ = (26) 



If we write yi^ = sin^ ^S—sin^ |-/3, y^ vanishes initially and the equation for y^ may be 
written 

4.y';' = Q' {h;' (cos ^-y;') -y^- h, sin /3 

+ (sin^i^+^i^)(cos'^i/3-^i^)[l/s-4g(sin^i^ + ?//)]}. . , . (27) 

We identify (27) with the equation 

iy';' = AqQ'{h,' + y;'){a,'-y;'){f;'-yn, (28) 

in which a^ = siv? ^ol—biv? \p. Equations (10)-(13) retain their form, and (7)-(9) 
become 

{by sin J^ + cos -J/3)' = 4g {K' + cos' i/3) (// - cos' J/S) cos' |a, . , . (29) 

6' = [h cos i/S-sin J/3)' = 4g (V-sln' 1/3) (// + sin' ^0) sin' Ja, . . . (30) 
-?>i'+(l/s) cos /3-l-4g sin' 1-^(2 cos' J/3 -sin' J-^) = ^q{a;'fi'-K'a^--htfi\ ■ (31) 

'•' T is the time interval between a zero or minimum and an adjacent maximum of the yaw. 
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The solution is 



l-g snu '' 



or 



in^i^^-cin^ i« - (sin'i«-sin'|-/3)cn^.t , . 

Ill qO — bill op — n n — ' "". . . . . . . \Ol^) 



sill 



The values of s and q^ in terms of A, sin |a, /c and g will differ from the values they 
had before by terms of the order sin^ \^ or & sin |/3 which are negligible in practice 
compared to a^}. Hence the previous solution may be applied provided that (32) 
replaces (10) in calculations of the curve of sin \^ as a function of Qi. For convenience 
in computing, (32) may be put in the form (neglecting g^) 

sin 1^ rr: ^^^ ^^ ^^^ X ^ SJU \^ SJU X ^ / (33). 

. ^ sin Q cos Q 

where 

tan = sin \a cot x/sin |-^, cos x = en t^. 

Formula (17) remains a valid approximation for h^ provided 6^ is defined by (30). 



§ 6. J^ Discussion of the Probable Effects of Damping, and Other Factors Omitted in the 

Foregoing Solution. 

Up to this point we have assumed that the motion in yaw is exactly periodic with 
half-period QT. This would be exactly true if the couple M were a function of S only, 
OP a fixed straight line, and no other couples existed. In actual fact, M is a function 
of the forward velocity and therefore of the time ; OP changes direction under the 
influence of gravity and the cross-wind force, and other couples besides M act on the 
shell, depending on the angular velocity of the shell. All these factors cause progressive 
changes in the curve of yaw from period to period ; for the case of the stable rounds 
they have been discussed at length in (A), where it is shown that they do not appreciably 
affect the determination of M at any velocity for small values of §. In particular, 
the effect of gravity is almost entirely allowed for by using (as we do) the true yaw 
and not the angle between the axis of the shell and some fixed straight line. As explained 
in (A)* the shape of the hole in the cards determines the true yaw and not the angle 
betAveen the axis of the shell and the normal to the card. 

There is no reason to expect that any of these damping effects will be relatively 
more important for an unstable than for a stable shell, except for the change of M with 
the velocity. Although the change in velocity over a single period is always small, 
yet when s is less than or nearly equal to unity a small change in M will cause a fairly 

* Log, cit., p. 318, footnote. 

Zi Ij £i 
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large change in the type of motion. This is shown clearly in fig. 3^ round III. 16, 
where the decrease in amplitnde and the change from an unstable to a stable type is 
quite marked, but still not sufficient to introduce any error in the determination of 
the couple for a single half -period. The effect is illustrated by the change in s in 
successive half -periods (Table III.), which is in general in the direction, and roughly of 
the amount, required by theory. 



0-45 




Kg, 3. The observed and calculated motion in yaw^ compared for selected rounds. The plotted 
points show the observed values of sin |S plotted against 0^ for rounds III. (11, 16) and IV. (8, 9). The 
continuous curves are the result of calculations described in detail in §§ 6, 8. Short vertical lines mark 
the positions of maxima and minima, and the origins of co-ordinates. The values of the constants used 
are as follows :— 

III. 11. (a) K = S0 degrees ; (6) k = 85 degrees, ^2^0. (c) k = 85 degrees, g^ ~ ■— 0*185. 

III. 16. (a) K = 80 degrees ; (b) k = 60 degrees, ^ = ; (c) k = 60 degrees, sin |/8 = 0-037 ; 
(d) K = 40 degrees, sin |/? == 0-037 ; (e) k = 0, with third order contact with (a) at maximum. 

IV. 8. (a) K =: 85 degrees ; (b) k ^ 75 degrees, sin |-^ === 0-034. 
IV. 9. (a) K ~ 60 degrees ; (5, c) k — 70 degrees. 

It appears that a change in M with v cannot alter an initial rosette motion into one 
with non-zero minimum yaw. This alteration, as in the stable case, must be due to 
the other couples depending on the angular velocity of the axis, and to the sideways 
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motion of the centre of gravity, whicli function as damping forces as explained in (A).*^ 
No means of dealing with these effects theoretically has yet been devised for the case 
of large yaw ; the observed changes are clearly of the general type which one's experience 
of the stable case would lead one to predict. 

§ 7. The Motion in cj>. 

It is not difl&cult to write down a formal expression for the motion in ^. If we take 
equation (25) and substitute for sin |^ from equation (32) we obtain after reduction 

, / _ Q , hQ sin 1^ r y 

^ ~ 1 + cosr (1-fcos^)^^' ^^ : 

— Q _, &1^ sin ^^ 1—g^m^u 

1 + cos (5 1 + cos 8 ( 1 —g^) sin^ |-/3 sn^ u + sin^ ^a cn^ t6 

where 6 is defined by (30). Now even when a is as big as 30 degrees there is still only 
a maximum difference of 7 per cent, between 1+cos 8 and 2, and this maximum is 
only effective for a short part of each period. Hence, for almost all purposes, we are 
still justified in replacing 1+cos (^ by 2t. In order to integrate this equation we notice 
that \Qdt = —du and that ^ == or t^ ::=: K corresponds to the minimum. Thus, 

. - . ■ in. . ^sin|-^ f^ jl-g'sn' u)du . ^ 



2X 



^ u 



{l—g^) sin^ ^/3 sn^ u + sin^ ^cc cn^ u 

This equation contains an elliptic integral of the third kind which can be evaluated 
in 0-f unctions. We have not, however, made this evaluation or calculated any actual 
^-curves from (34) mainly because it does not appear that any further information as 
to the forces acting on the shell would be obtained thereby. We shall content ourselves 
in this paper with a statement of sufficient theoretical results to show that the observed 
0-curves are qualitatively of the form to be expected from (34). A more detailed 
discussion of these curves, however, would, we think, be of some interest. 

It is convenient to treat the motion by using the variable i/^ ^^^ 95>— 95>o— 10^. 
When /3 is zero, ^' will be constant and equal to |1} to our present approximation ; 
with y and ^r as polar co-ordinates the motion then consists of an oscillation in a straight 
line through the origin, for i^' ^ o. In the general case we may eliminate dt between 
equations (33)' and (28), and on substituting for b from (30) obtain an equation for 
dy/d\fr in the form J 

4iJ ^Vsin'i/8 /\ sin^WV /^i'-sin^i/5/\ fUsin'^/Bj' ' ^ ' 
If we assume that hi^ and /^^ are large compared to sin^ |a, the last two brackets 

'''•' In particular see p. 313. 

"f It is easy to estimate the precise efiect of this approximation in the simple case of the rosette motion. 
The error caused is always very small. 

J In deducing (35) we replace 1+cos 8 by 2. 
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in (35) are effectively unity, and the equation reduces to that of an ellipse in polar 
co-ordinates with axes sin |a, sin |-/3. This tends to a straight line as the limiting 
form when /3 ->0. 

The same (simplified) relation between § and yjr was obtained in (A)* for the case 
of small yaw only. The elliptic motion may be calculated most easily by recognising 
that xfr is the auxiliary angle of equation (33) so that 

<P = (po + i^^ + 0. 

The conditions under which this approximation is valid may be seen by reference to 
equations (11) and (12) which are satisfied hj fx, lfi\ and a^. We notice that/^,^ is 
never small compared to unity and tends to infinity as q^ and therefore g tends to zero ; 
while A^^ may be comparable with a^ unless h is small. Thus the approximation is 
really only valid for <§ > 1, in which case it applies even if ^ is not small compared 
to a. Finally, if 1— A^^ is small, equation (12) shows that li^ will be small compared 
to a^^ and the third bracket in (36) will be more important than the second near a 
minimum of y ; this indicates that the shape of the {y, \f/) curve there approximates 
to an hyperbola instead of to an elongated ellipse ; the curve may also no longer be 
re-entrant, the total change of \fr in one period differing from ±7r. 

Examples of all these results may be seen in fig. 4, in which the observations for 
three different rounds are plotted with {y, \/y) as polar co-ordinates. For round I. (5), 
for which /3 is small and the shell just stable, we find the expected elongated ellipse-like 
curve, with a slowly-developing minor axis caused by the damping factors. In III. (16) 
/c = 80 degrees for the first half -period, diminishing to 40 degrees for the third ; a 
considerable minimum yaw developes and the curve is less like an ellipse, though the 
maxima are still nearly 180 degrees apart. Finally, in IV. (8) k == 85 degrees, falling 
to 75 degrees, and the shape of the curve near the minimum clearly resembles an 
hyperbola ; we may guess then the angle between the two maxima is somewhere about 
100 degrees instead of 180 degrees. 

Lastly, a word must be said about the observational determination of Q, which is, 
of course, theoretically determined by the muzzle velocity, twist of rifling and moments 
of inertia of the shell. In all cases the slope of the ^-curve over the first half-period, 
or rather more, is uniform and well determined. Since /3 appears to be really very 
small initially one may expect from theory the slope of this part of the curve to be 
|0 whatever the value of k. The agreement between this observed slope and the 
calculated value of Q is satisfactory. 

§ 8. The Method of Analysis. 

The method of analysis of the sin |(5-curves will now be explained with reference 
to fig. 3. After the observed values of sin ^§ have been plotted against Q^, the values 

* Log. cit,, p. 346, equation 4.06. 
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of Qt corresponding to the ends of half -periods can be determined from the symmetry 
of the curve, with the exception of the start of the first half-period. This must, of 
course, be near the muzzle, but need not be actually at the muzzle, and its precise 
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Fig. 4. Curves showing the observed angular motion of the axis with (sin JS, \p) as polar 
co-ordinates {\p = (jy~4>Q—^^t). The circles show the observed points and the curves are 
drawn freehand through them. 

For I. (5) 1 cm. represents 0-01 in sin|S. 

For III. (16) and IV. (8) 1 cm. represents 0-02 in sin J8. 

position must be guessed. The values of sin ^a can be obtained from a rough curve 
drawn freehand through the observations. Each half-period is then analysed separately 
and the only additional constant (when ^ — 0) required for computing a curve, by 
the approximate formula 

sin ^§ = sin |-a en (K — XQt, k)^ 

is the value of h or sin k. This may be approximated to with the help of the following 
artifice. Draw by eye a cosine curve (fig. 3, III. 16, first half-period) to have third 
order contact with the observed sin i^ (5- curve at the maximum and to cut the axis 
at T. Then the ratio of NT to the half -period NP must be 7tj2K.. For all curves of 
the form i?/ — A en {K.—j3x), for different values of h, have third-order contact at the 

A first 



nUT> m 



maximum ordinate, and a cosine curve is the limiting case, in which K - 

estimate of k can usually be made by this method to the nearest 5 degrees, when 
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AC > 45 degrees. If necessary, similar curves are drawn for other values of k (rounds 
III. 11, III. 16) and the true values of k, sin|a and OT (the length of the half-period) 
finally settled by interpolation. 

If required a curve can be calculated by the exact formula (10) after g has been 
determined by (19), but the change of shape is negligible unless sin |a and, therefore, g 
is very large. An example of the effect of including g is shown in fig. 3 for III. (11), 
but round III. (11) and its companions are the only ones in which g had any sensible 
effect. Rounds IV. (8 and 9) illustrate the effect of a considerable alteration in b 
(representing the initial disturbance) in causing, when s < 1, a considerable change in 
period but only a small change in amplitude. The fit obtained between calculated 
and observed curves is generally good. The selected curves in fig. 3 are a fair sample 
of the whole. 

When the minimum yaw j3 is not zero the curve is first calculated as above, as if 
^ =^ 0, and then corrected by (33), using the observed value of /3, which is obtained 
like a from a rough curve. An example of such curves will be found in fig. 3 for the 
second half-period of IV. (8) and the second and third half-periods of III. (16). After 
the 13 correction has been put on, the value of k may require readjustment to obtain 
the proper fit. 

The values of k, sin |a and QT so determined are given for each half-period of each 
round in Table II., together with sin|/3 and X obtained from the equation \ ~ K/iiT. 
We then obtain s, q and h from equations (15)-(17).'^ From the values of s and q and 
the other observational data we can calculate /^ sin § as a function of v/a and S from 
formulae (3) and (4). The results are shown in figs. 1 and 2, and Tables I. and III., 
and have already been discussed. 

The damping effects appear in the variation of the various constants from one 
half -period to another. In general, s increases approximately at the rate required by 
theory, i.e., inversely as the square of the velocity at the middle point of the half- 
period." 



* When K is mucli less than 45 degrees the method breaks down, as k cannot be determined satisfactorily 
from the observational curves. The method explained in (A), p. 343, could then be employed. This is 
equivalent to assuming ^ = and using formula (16) to determine k given sin Ja and A. Under these 
conditions sin Ja is so small that the value of q does not appreciably aifect either the value of s or the shape 
of the curve of /m against 5. 

f Theoretically, s should increase while q should remain roughly constant. But sin |a is common to 
the first and second half -period, while A is determined mainly by the shape of the curve near the maximum. 
Hence, in general, k alone varies between the first and second half -periods. It appears that the result of 
varying k only in formulae (15, 16) is to produce a fictitious decrease in q while the increase of s is diminished 
as may be seen in Table II. For this reason mean values for the whole period are used in constructing 
figs. 1 and 2. 
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Table II. — General Table of Eesults. 



Column 1. Number of round and number of half -period for the round (in brackets). 

2. Muzzle velocity, f.s., for round or mean for group. 

3. Air density p, lb./(ft.)^, and temperature^" F. 

4. Q {=^ AN/B), radians/sec. 

5. QT, radians, where T is the observed duration (sec.) of each half-period. 

6. Observed values of sin |a, where a is the maximum yaw. 

7. Observed values of sin J/3, where /3 is the minimum yaw. 

8. Values of /c, degrees, where y^ (~ sin /c) is the modulus of the elliptic 
function which fits the observations. 

9. Values of x(" K/i^T), where K is the complete elliptic integral of the 
first kind for the value of k in (8). 



J) 



>5 



53 



?3 



>> 



?? 



:>:i 



1 


2 


3 


4 


5 

17-5 
19-8 
16-2 
15-2 


6 


7 


8 


9 


1.8 (1) 
I. 9 (1) 

(2) 
I. 10 (1) 


1049 
1084 

1084 


0-0792 
43° 


84-2 
86-9 

86-9 


0-276 
0-251 
0-251 
0-245 




79 

83 
70 

75 


0-172 
0-154 
0-155 

0-182 

i 


1.17 (1) 

(2) 
I. 18 (1) 

(2) 


1312 


0-0812 
40° 


105-2 


16-6 
13-8 
13-5 
12-0 


0-192 
0-192 
0-203 
0-203 


0-020 
0-049 


75 
65 
75 
70 


0-167 
0-167 
0-205 
0-209 


1.5 (1) 

(2) 

(3) 
I. 6 (1) 

1 (2) 

I. ^ (1) 

(2) 


1582 


0-0782 
45° 


125-4 


13-8 
15-6 
12-0 
13-6 
14-4 
24-0 
19-5 


0-108 
0-108 
0-070 
0-118 
0-118 
0-090 
0-090 


0-018 


40 

57 
20? 
45 
35 
56 
0? 


0-130 
0-133 
0-135 
0-136 
0-120 
0-086 
0-082 


11.11 (1) 

(2) 
II. 12 (1) 

(2) 
II. 13 (1) 

(2) 


1107 


0-0807 

42° 


81-9 


18-2 
14-6 
20-8 
16-4 
22-0 
14-4 


0-161 
0-161 
0-148 
0-148 
0-183 
0-183 


0-055 
0-075 
0-046 


57 

20? 

52 

35 

70 

53 


0-114 
0-112 
0-095 
0-105 
0-114 
0-138 




1 

i 

i 
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Table II. (continued) 



1 


2 


3 


4 


5 


6 


7 


8 


9 


11. 14 (1) 


1334 


0-0812 


98-6 


18-4 


i 

0-135 




60 


1 

i 

0-119 1 


(2) 




40^ 




19-8 


0-135 


0-014 


55 


0-105 1 


IL15\(1) 








16-6 


0-122 




48 


0-117 1 


* 16/(2) 

1 








15-9 


0-122 


0-012 


35 


0-111 1 

! 
1 

i 

1 


III. 11(1) 


1077 


0-0807 


99-0 


34-0 


0-332 


j — 


1 

84 


■ ■ ' i 

i 

0-107 


111.12(1) 




42° 




32-0 


- 290 


1 


85 


0-120 


111.13(1) 








28-4 


0-301 

i 

[ 


i 


85 


0-135 

i 
! 

f ' 
1 I 


111.14(1) 


1312 


0-0812 


120-4 


22-0 


0-204 




78 


i 

0-135 


(2) 




40° 




17-2 


0-204 




70 


0-146 


III. 15 (1) 








18-0 


0-207 




75 


0-154 


(2) 








14-5 


0-207 


0-053 


55 - 


0-141 


(3) 








13-9 


0-187 


0-053 


55 


0-147 


111.16(1) 








21-8 


0-204 




80 


0-145 


(2) 








15-6 


0-204 


0-037 • 


60 


0-138 


(3) 








17-0 


0-179 1 


0-037 


40 


0-105 


IV. 10 (1) 


884 


0-0807 


62-5 


16-5 


0-233 ; 


) 


70 


0-152 


(2) 




42° 




13-3 


0-233 : 


— , 


50 


0-145 1 


IV. 11 (1) 








21-0 


0-225 : 




78 


0-142 


(2) 








15-0 


0-225 ^ 


— 


50 


0-129 


IV. 12 (1) 








20-0 


0-207 : 


— 


77 


0-145 


(2) 








14-5 


0-207 1 


i 


62 


0-153 


IV. 7 (1) 


1553 


0-0779 


109-7 


1 

1 
( 

15-5 


0-213 J 


i 


i 
83 


i 

0-226 


IV. 8 (1) 




47° 




18-0 


0-199 





85 


0-213 


(2) 








13-0 


0-199 


0-034 


75 


0-213 


IV. 9 (1) 








11-5 

1 


0-218 




70 

! 


0-218 


IV. 1 (1) 


2130 


0-0782 


150-6 


12-5 


i 

0-160 


~ 


66 


1 

0-187 


(2) 




45° 




12-9 


0-160 


0-018 


65 


0-179 1 


IV, 2 (1) 








18-0 


0-157 




82 


0-187 


(2) 








14-2 


0-157 


0-013 


70 


0-176 1 


IV. 3 (1) 








17-0 


0-160 




80 


0-185 


(2) 








14-8 


0-160 


— - 


75 


0-187 


IV. 5 (1) 








15-5 


0-160 




74 


0-175 ! 


(2) 








13-9 

i 
! 


0-160 


0-015 


70 


0-180 



* II. (15 and 16) are practically indistinguishable. These values refer to their meah. 



THE AERODYNAMICS OP A SPINNING SHELL. 



245 



Table III. — Mean values, for each group, of s, q and v ; the corresponding values of 
v/a and /m {vja, 0) ; and the percentage spread in s. Each half -period is treated 
separately in each group. 



Group. 


s. 


Percentage 

spread 

in s. 


^- 


V. 


vja. 


Jm- 


I. 8-10 (1) 
L9 (2) 
I. 17, 18 (1) 
I. 17, 18 (2) 
I. 5-7 (1) 
I. 5 7 (2) 
I. 5 (3) 


0-896 
0-932 
0-892 
0-908 
1-000 
1-007 
1 - 060 


0-8 

4-3 

5-5 
2-3 
5-6 


0-213 

0-084 
0-578 
0-597 
0-398 
0-507 


1063 
1053 
1297 
1264 
1560 
1517 
1494 


0-965 
0-956 
1-181 
1-150 
1-413 
1-375 
1 - 354 


11-80 
11-56 
11-61 
12-01 
10-56 
11-08 
10-86 


II. 11-13 (1) 
II. 11-13 (2) 
II. 14-16 (1) 
II. 14 16 (2) 


0-977 
1-014 
0-987 
1-006 


3-1 

5-5 
2-1 
3-2 


0-068 
0-04 
0-284 
0-066 


1094 
1075 
1315 
1286 


0-994 
0-977 
1-198 
1-171 


10-34 
10-30 
10-20 
10-47 


III. 11-13 (1) 
III. 14-16 (1) 
III. 14 16 (2) 
III. 15, 16 (3) 


0-946 
0-931 
0-959 
-,990 


2-6 
1-8 
3-6 
3-7 


0-119 
0-203 
0-118 


1063 
1296 
1263 
1239 


0-966 
1-180 
1 - 150 
1-128 


13-08 
13-10 
13-43 
13-52 


IV. 10-12 (1) 
IV. 10-12 (2) 
IV. 7-9 (1) 
IV. 8 (2) 
IV. 1-3, 5 (1) 
IV. 1-3, 5 (2) 


0-934 
0-975 

0-852 
0-864 
0-897 
0-907 


0-3 
3-8 
4-5 

3-5 
3-5 


0-146 
0-033 
0-793 

0-820 
0-987 
0-885 


877 
866 
1532 
1488 
2102 
2052 


0-797 
0-787 
1-385 
1-345 
1-905 
1 - 860 


10-00 
9-79 
11-40 
11-93 
10-80 
11-20 



The values of /m here given are not corrected for the effect of the cards. 



§ 9. The Values of b. 

The value of b represents the initial angular velocity of the axis of the shell, for at the 
beginning of the first half-period rosette motion may be assumed and §' = &Q. The 
values of b for each round are given in Table IV. 

Table IV, — Values of b { = S'olQ), where S^^ is the initial angular velocity of the axis 

of the shell. 



Group. 


Values of h. 


Mean. 


Group. 


Values of b. 


Mean. 


I. 8-10 


0-018 


0-011 


0-023 


0-017 


III. 11-13* 


0-007 


0-006 


0-007 


0-007 


I. 17, 18 


0-017 


0-022 




0-019 


III. 14-16 


0-011 


0-016 


0-011 


0-013 


1.5-7 


0-021 


0-023 


0-009 


0-018 


IV. 10-12 


0-024 


0-013 


0-013 


0-017 


IL 11-13 


0-020 


0-017 


0-014 


0-017 


IV. 7-9 


0-012 


0-007 


0-032 


0-017 


II. 14-16 


0-016 


0-018 


0-018 


0-017 


IV. 1-3, 5 


0-024 


0-008 


0-010 
















0-015 






0-014 



2 M 2 
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They represent the size of the initial disturbance which upsets the nose-on motion, 
and vary irregularly in any one group, as might be expected. Their general size, 
however, is remarkably consistent from group to group. The corresponding values 
for the stable rounds and rough values for some of the rou.nds here analysed have been 
discussed by us in a previous paper.* The further results here given confirm the state- 
ments of that paper ; the mean value of b for any group is practically constant some- 
where between 0-01 and 0-02 for all groups fired. 

§ 10. Alloivance for the Effect of Cards, 

In obtaining figs. 1 and 2 and Table I., but not in other cases, a small correction 
has been made to the value of the couple to allow for the impulsive action on the shell 
when it strikes a card. The amount of this correction was calculated from a few special 
rounds fired with cards on the far screens only. The following argument indicates 
that the effect should be roughly proportional to the couple due to the air and 
independent of other factors except the number of cards and the muzzle velocity. 
The effect of each card may be considered as an impulsive force whose moment about 
the centre of gravity is roughly proportional to sin ^. If the cards are (as they were) 
uniformly distributed in space and therefore in time, there will be a total impulse per 
second which is proportional to sin §^ and if the cards are not too far apart this is 
equivalent to a steady couple roughly proportional to the couple due to the air. The 
correction worked out in the case of the stable rounds at from 3 per cent, to 4 per cent, 
depending on the muzzle velocity. In view of the preceding argument the same 
corrections are applied here. 

§ 11. Concluding Remarks, 

In view of possible future experiments it is of interest to compare the merits of the 
stable and unstable rounds for the purpose of this analysis. A large proportion of 
the rounds analysed in (A) had a stability coefficient s of about 1 • 8. The advantages 
of this are that the theory of the motion is practically complete, and in addition to 
the values of f^^ as a function of vja rough values of the damping factors were obtained, 
which could be greatly improved if a longer range were available. The maximum 
yaw, however, is small, so that small errors in the determination of the yaw are 
important, while the periods are short so that, unless the cards are close together, it is 
difficult to draw curves through the observation points with sufficient certainty to 
determine the periods accurately. 

For the unstable rounds here analysed the theory is imperfect as regards the 
determination of damping. But a small change in the value of the couple M will 

* ^Proc. Camb. Phil. Soc./ vol. XX, p. 31L 192L 
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make a large change in the type of motion, so that the method is now a very sensitive 
one for determining the conple, and no great accuracy of observation is required. 
The yaw developed is also much larger, so that values of the couple are determined 
over a larger range of yaws. 

From the point of view of general aerodynamical theory the results form a preliminary 
contribution to the problem of determining the force system impressed by the air on a 
body moving unsymmetrically through a fluid at velocities at which the compressibility 
of the fluid produces marked effects. 



